The temperature-dependent elastic properties of B2-MgRE 
intermetallic compounds from first principles * 



Rui WangJ Shaofeng Wang, and Yin Yao 
Institute for Structure and Function and department of physics, Chongqing University, 
Chongqing 400044, People's Republic of China. 

Abstract 

Using the density functional theory (DFT) formulated within the framework of the plane-wave basis pro- 
jector augmented wave (PAW) method, the temperature-dependent elastic properties of MgRE (RE=Y, Dy, 
Pr, Sc, Tb) intermetallics with B2-type structure are presented from first-principles. Our calculations are 
based on the fact that the elastic moduli as a function of temperature mainly results from thermal expansion. 
The comparison between the predicted results and the available experimental data for a benchmark material 
NiAl provides good agreements. At T — OK, our calculated values of lattice parameter and elastic moduli 
for MgRE intermetallics show excellent agreement with previous theoretical results and experimental data. 
While temperature increases, we find that the elastic constants decrease and approach linearity at higher 
temperature and zero slope around zero temperature. 

PACS: 71.20.Lp, 62.20.D-, 65.40.De, 71.15.Mb. 
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1 Introduction 

Since magnesium (Mg) alloy often shows high strength to weight ratio JT), it is one of the potential materials for 

technological applications. Especially, Mg alloys have been attractive for the applications of aeronautical and automotive 

industry. However, the low strength and creep resistance at high temperature is a very serious problem in application of 

magnesium alloys. Therefore, it is quite necessary to improve the strength and creeping resistance of Mg-based alloys 

at elevated temperatures. Recently, it has been reported that some rare-earth-Mg intermetallic compounds MgRE (where 

RE indicates a rare-earth element) with B2-type structure have good creep and high temperature strength [HOB). MgRE 
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intermetallics are extremely attractive structural materials for applications in automobile parts and aerospace industries. 
So various properties such as magnetic properties, linear and nonlinear elasticity, stacking fault, and thermal properties 
for the B2-MgRE intermetallics have been intensely investigated iE] [6l d [8] |9] [TOl [TT1 |T2l |T3l [141 . In our previous study, 
the temperature dependence of various quantities such as the thermal expansions, bulk modulus, and the heat capacity are 
obtained by using the density functional theory (DFT) and density functional perturbation theory (DFPT) IfTTI . 

In general, elastic properties of a solid are very important because they are closely associated with various fun- 
damental solid-state properties such as interatomic potentials, equation of state, and phonon spectra. The temperature 
dependence of the elastic constants of a material is important for predicting and understanding the mechanical strength, 
stability, and phase transitions of a material [ 15 1. Though modern electronic-structure methods based on DFT can treat the 
zero-temperature pressure dependence elastic moduli from first principles accurately lfl6l . treating the corresponding tem- 
perature dependence of elastic constants is still a formidable challenge IfTTI . However, the temperature-dependent elastic 
properties have also been investigated by using the first-principles DFT methods in several groups |[T5lfl^fT8l[T9ll20ll2TI . 
Especially, Wang et. al. Il20l |2T1 demonstrate that the temperature dependence of elastic moduli mainly results from 
volume change as a function of temperature, and their applications to many materials show excellent agreement between 
the calculated values and experimental data. In this paper, we employ the first-principles method 11201 |2T1 to study the 
temperature-dependent elastic constants for MgRE (RE=Y, Dy, Pr, Sc, Tb) intermetallics with B2-type structure. To judge 
that our computational accuracy is reasonable, we have also computed the temperature dependence of elastic moduli for 
a benchmark material NiAl, which has the same B2-type structure and accompanied by available experimental data. 

2 Theoretical methods 

The Helmholtz free energy of intermetallics at a constant volume V and T has three major contributions Il22l|23ll 



where fistic is the zero-temperature energy of a static lattice, F e \ is the thermal electronic contribution to free energy from 
finite temperature, and F p h represents the phonon free energy arising from the lattice vibrations. We obtain both fistic an d 
F e j from first-principles DFT calculations directly, assuming that the eigenvalues for given lattice and ion positions are 
temperature independent and only the occupation numbers change with temperature through the Fermi-Dirac distribution 
lfT9ll24l . The phonon free energy within the framework of quasiharmonic approach (QHA) is obtained from 



F{V, T) = £ static ( V) + F d ( V, T) + F ph (V, T), 



(1) 




(2) 



where a> K represents an individual phonon frequency. 
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The isothermal elastic constants are derived from the second-order strain derivative of the Helmhotz free energy [25 

1 d 2 F 



(3) 



T.E' 
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where s' indicates that all other stains are held fixed. The elastic behaviour of B2-MgRE intermetallics is completely 
described by three independent constants C\ v C T n , and Cj 4 (in Voigt notation). The bulk modulus is defined by a linear 
combination of the elastic constants 

B T = (Cf, + 2C[ 2 )/3, (4) 



and determined from the Vinet equation of state [26 1 corresponding to a pure volume deformation of the lattice. We apply 
a tetragonal strain to calculate {C\ x - C T l2 ) 



(60 
6 
(1+6T 2 - 1 ) 



(5) 



and the corresponding free energy is 



F(V, 6) = F(V, 0) + 3(Cfi - C T l2 )V6 2 + 0{6\ 
where F(V, 0) is the free energy of the unstrained state. Cj 4 is obtained from a volume-conserving monoclinic stain 



(6) 
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which leads to the free energy change 



F(V,S) = F(V,0) + 2C T AA V6 2 + 0(6 4 ), 



(8) 



We calculate the free energy for different strains, and then the elastic constants are calculated from the fitted quadratic 
coefficients. To compare with experiment, the isothermal elastic moduli must also be transformed to the adiabatic elastic 
moduli by the following relation E71 



^ v kk> 



(9) 



where CV is the specific heat at constant volume and is the linear thermal expansion tensor. For B2-MgRE intermetallics 
with cubic crystal, Eq. Q simplifies to IfTTll 



C s = C T 

*-44 "-44' 



(10) 



with a is the thermal expansion coefficient. From Eqs. ([4]) and ( 10 1, it follows that the adiabatic bulk modulus is just 

TV , , 
B S =B T + —a 2 B\ 

Cv 



(11) 



We use a quasistatic approach, which was recently developed by Wang et al J20l I2H . to calculate the temperature- 
dependent elastic constants Cij(T). In the quasistatic method, it is assumed that the temperature dependence of elastic 
constants is solely caused by thermal expansion. The applications of quasistatic approach to elastic constants of many 
materials from OK up to their melting points show excellent agreement between the computed values and experimental data 
1 20 21 1. In this paper, the temperature-dependent elastic constants are obtained in the following three steps [20|. The first 
step is calculating the thermal expansion and the equilibrium volume V(T) at T by using the first-principles quasiharmonic 
approximation, while isothermal bulk modulus Bj as a function of temperature T is also obtained from fitting to the Vinet 
equation of state l26l . The detailed procedure of calculating the thermal expansion coefficient, isothermal bulk modulus 
Bj were shown in our recent work [ 1 1 1. In the second step, we predict the static elastic constants as a function of volume 
Cjj(V) by using the energy-strain relation based on Eqs.(|4jj8j. In the third step, the calculated elastic constants from the 
second step at the volume V(T) are approximated as those at finite temperatures, i.e., Cf^T) = Cf.[T(V)]. Then, the 
adiabatic elastic constants Cfj can be obtained from Eqs. 

The computational approach is based on the density-functional theory (DFT) and density-functional perturbation 
theory (DFPT) as implemented in the VASP package j28l [29] EE |32 • We use the Perdew-Burke-Ernzerhof (PBE) |32] 
33 1 generalized gradient approximation (GGA) for the exchange-correlation functional. The plane-wave basis projector 
augmented wave (PAW) method l34l [35l is used. Since high accuracy is needed to evaluate the elastic constants, the 
convergence of strain energies with respect to the Brillouin zone integration was carefully checked by repeating the 
calculations for 21 x 21 x 21 and 25 x 25 x 25 Monkhorst-Pack [36 1 meshes, and we found that fluctuation both for C44 
and n is lower than 0.5GPa. Hence, we used 21 x21 x21 in the full Brillouin zone giving 726 irreducible k-points. In 
addition, we used a high plane-wave energy cutoff of 600eV which is sufficient to calculate the elastic moduli accurately. 

To calculate the phonon frequency, we have carried out DFPT calculations on this 3x3x3 supercell with 54 atoms 
by using PBE-GGA exchange-correlations effects and 7x7x7 k-point grid meshes for Brillouin zone integrations. 
The phonon free energy were obtained by using PHONOPY 1371 [38l |39ll package which can support VASP interface to 
calculate force constants matrix directly. To consider the thermal expansion, we calculate the vibrational frequencies 
of all calculated interme tallies at 13 volumes. The thermal expansion coefficients a and the thermal bulk modulus B T 
are obtained by minimizing Helmholtz free energy Eq. ([T} with respect to V from fitting the integral form of the Vinet 
equation of state (EOS) |26|. 

3 Results and discussion 

The calculated equilibrium lattice constants a, elastic constants C,y, and bulk modulus B at T — OK for MgRE (RE= Y, 
Dy, Pr, Sc, Tb) in comparison with the previous calculated results J5] [Toll are listed in Table [T] Where available, the 
experimental data [40] for lattice constants are also listed for comparison. The calculated lattice constants are in excellent 
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agreement with experimental data within 0.5%. While the lattice constants in present calculations shows no difference 
from the results obtained by Tao et. al. lflOl . since we employ the PBE-GGA exchange-correlation functional compared 
to PW91-GGA in their calculation. No temperature-dependent elastic constants have been reported theoretically and 
experimentally for all MgRE intermetallics. At T — OK, our calculated results for elastic constants agree well with the 
previous DFT data (9l [Toll and it is noticeable that, the requirement of mechanical stability for all calculated MgRE, 
namely, C u - Cu > 0, C u > 0, and C44 > (42), is satisfied. 

The calculated isothermal and adiabatic bulk modulus for NiAl and MgRE (RE=Y, Dy, Pr, Sc, Tb) as a function 
of temperature are shown in Figure [T] The all isothermal bulk modulus Bj are obtained by minimizing Helmholtz free 
energy Eq. ([T]) with respect to V from fitting the integral form of the Vinet equation of state (EOS) |26|. Figure[Tfa) shows 
our findings for benchmark materials NiAl, which has the same B2-type structure as MgRE intermetallics, accompanied 
by available experimental data form Rusovic and Warlimont ATI , The calculated values of adiabatic bulk modulus B$ 
for NiAl are in good agreement with experiments. Our results also show no discrepancy with previous DFT calculated 
values obtained from using the quasistatic approach by Wang et. al. GUI . The overall observation is that both the Bj 
and Bs decrease with increasing temperature and approach linearity at higher temperature and zero slope around zero 
temperature. It it worth to note that the corrections TVa 2 Bj/Cy between adiabatic elastic moduli and isothermal ones 
increase with increasing temperature. Among the five intermetallics [shown in Figure [TJb)-(f)], through the temperature 
range, MgSc [Figure [TJe)] and MgPr [Figure [TJd)] have the highest and the lowest bulk moduli, respectively, i.e., MgSc 
is the most incompressible and MgPr is the most compressible. 

Figure [2] shows that the calculated temperature dependences of the adiabatic elastic constants , C s n and C^ 4 as a 
function of temperature. In Figure |2j (a), we compare the calculated values of C^. for NiAl with experiments BTl! and get 
a useful test of the accuracy of the method and the precision of our calculations. For all calculated MgRE intermetallics, 
we find the elastic constants Cf. decrease with increasing temperature, since thermal expansion may soften the elastic 
moduli at high T. We also find that the trend of Cf. approach linearity at higher temperature and zero slop around zero 
temperature. The calculated curves from Figure [2] indicate that the elastic constants of MgRE follow a normal behavior 
with temperature. For MgY [Figure|2^b)], MgPr [Figure[2jd)], and MgSc [Figure^e)], we find that the curves of , C s n 
and are nearly parallel, i.e., almost decrease to same extent in the whole temperature range. For MgDy [Figure |2jc)] 
andMgTbpiguregf)], C s n decrease to slightly larger than C s n and C^. This character indicates that MgRE intermetallics 
can keep their mechanical properties through wide temperature range and have good high temperature stability. 

4 Conclusions 

In summary, we have performed first-principles quasistatic approach to study the temperature-dependence of the elas- 
tic moduli of five MgRE (RE= Y, Dy, Pr, Sc, Tb) intermetallic compounds with B2-type structure. Our calculations are 
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based on the fact that the elastic moduli as a function of temperature mainly results from thermal expansion. The elastic 
constants as function of volume are determined by the first-principles DFT total-energy calculations within the framework 
of the method of homogeneous deformation and the thermal expansion is determined by first-principles phonon calcula- 
tions based on DFPT with QHA. To benchmark the reliability results of the presented method, the comparison between 
the predicted results and the available experimental data for a benchmark material NiAl provides good agreements. At 
T - OK, our calculated values of lattice parameter and elastic moduli for MgRE intermetallics show excellent agreement 
with previous theoretical results and experimental data. We find that the elastic constants follow a normal behavior with 
temperature: decrease with increasing temperature and approach linearity at higher temperature and zero slope around 
zero temperature. 
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Table 1 : Our calculated lattice constants a, elastic constants Cy, and bulk modulus B for B2-MgRE (RE=Y, Dy, Pr, Sc, Tb) at T = OK 
compared to previous computed results and experimental data. Note that Cf, = Cj. = Cy at T = OK. 





MgY 


MgDy 


MgPr 


MgSc 


MgTb 


a (A) 


3.795°, 3.7%", 3.795 1 ', 3.796" 


3.778", 3.765", 3.776°, 3.759° 


3.910", 3.901", 3.909°, 3.912° 


3.580°. 3.593°, 3.580°, 3.597" 


3.789", 3.781°, 3.789° 3.781° 


Cn (GPS) 


52.60", 53.37'', 53.07 r 


51.29", 51.86°, 51.85° 


50.61°, 49.96°, 49.45° 


64.93", 70.77°, 63.86° 


52.59", 53.32°, 52.46° 


Cn (GPa) 


35.58°, 36.39°, 36.10° 


36.38", 37.57°, 36.85° 


30.03°, 30.80°, 31.08° 


43.51", 43.11°, 45.09° 


35.58°, 36.18°, 36.20° 


C44 (GPa) 


38.12°, 39.05°, 39.26° 


38.63", 38.35", 39.74° 


37.10", 36.80°, 36.20 


53.61", 55.65°, 52.62° 


39.21°, 39.82°, 39.95° 


B (GPa) 


41.25°, 42.06°, 41.75° 


41.35", 42.32", 41.85° 


36.86", 37.20°, 37.20° 


50.65", 52.33°, 51.35° 


41.25°, 41.89°, 41.62° 



" This work; 

Ref. 9 from first-principles calculations; 
c Ref. 1 1 Q] from first-principles calculations; 
d Ref. 1 40| from experiment. 
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Figure 1 : (Color on line) The temperature-dependent adiabatic bulk modulus B s (solid line) as well as the calculated isothermal bulk 
modulus B T (dashed line) for (a) NiAl, (b) MgY, (c) MgDy, (d) MgPr, (e) MgSc, and (f) MgTb. For NiAl, the full symbols denote the 
corresponding experimental data form Rusovic and Warlimont I4TI. 



9 



cd 
Q_ 

a 200 

CO :=- 

o 

I 150 



1 100 

CD 

o 



CO 
X! 

< 



50 



• 


• 


(a) NiAl 
• 


• Sf 




■ 


■ 




pS 
°12 


• 


_ -♦- _ 
















44 















55 



200 



400 600 
Temperature [K] 



800 



1000 



CO 

CL 



a so 



i= 45 



CO 

o 

c 
ro 
to 

o 
o 
o 
to 
iS 

CD 



-O 

< 



40 



35 



25 



(b) MgY 



44 



12 



200 400 600 800 1000 
Temperature [K] 



55 
o 50 



o 



45 



CO 

| 40 

o 
o 

CO 

J5 

CD 

o 

CO 

I 25 
< 

20 



30 





(c) MgDy 








^ir 






pS 
°44 






"- c s- - -'-'=: 

12 









200 400 600 800 1000 
Temperature [K] 



CO 

CL 

O 



55 



50 



CO :=■ 

° 45 

CO HJ 

c 

CO 

I 40 

o 
o 

M 35 

CO 

iS 

CD 

o 30 
to 

JD 

55 

< 



25 



20 



'(d) MgPr 



12 



200 400 600 800 1000 
Temperature [K] 



70 
I 65 
co.. 60 
| 55 

CO 

1 50 
o 



45 
to 

1 40 



to 35 

CO 
< 



30h 
25 





.(e) MgSc 




CO 

Q. 






......... 


CD 






CO :=- 

o 

CO 














r s 


consta 








o 






C s 

_ _12 


ic elast 








iabat 








"O 

< 



200 400 600 800 1000 
Temperature [K] 




200 400 600 800 1000 
Temperature [K] 



Figure 2: (Color on line) The temperature-dependent adiabatic elastic constants for (a) NiAl, (b) MgY, (c) MgDy, (d) MgPr, (e) 
MgSc, and (f) MgTb. The solid, dashed-dotted, and dashed curves represents the values of C s lr C s l2 and Cf,, respectively. For NiAl, 
the full symbols denote the corresponding experimental values form Rusovic and Warlimont 1411 . 
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